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METHOD FOR RENDERING DIGITAL RADIOGRAPHIC IMAGES FOR 
DISPLAY BASED ON INDEPENDENT CONTROL OF FUNDAMENTAL 
IMAGE QUALITY PARAMETERS 

5 FIELD OF THE INVENTION 

This invention relates in general to image enhancement in digital 
image processing and in particular, to a method of image enhancement to be used 
in a medical radiographic imaging system, such as a computed radiography 
system. 

1 0 BACKGROUND OF THE INVENTION 

It is a common desire to enhance images acquired from imaging 
devices to improve their image quality. There are fundamental attributes that 
govern the image quality of a grayscale image. These attributes are the brightness 
of an image, the dynamic range of an image, the contrast of detail in an image (or 

15 detail contrast), the sharpness of edges in an image, and the appearance of noise in 
an image. It is beneficial to provide a system that enables direct and independent 
control of these attributes of image quality. Further a system, which enables the 
control of the detail contrast, sharpness, and noise appearance in a density 
dependent fashion is also beneficial. 

20 U.S. Patent No. 5,805,721, inventors P. Vuylsteke and E. 

Schoeters, entitled "Method and Apparatus for Contrast Enhancement", issued 
Sept. 8, 1998 describes a multi-resolution method for enhancing contrast with 
increased sharpness that includes dynamic range compression, and enhancing 
contrast without remarkably boosting the noise component. The described 

25 invention enhances detail contrast and sharpness via a multi-resolution method 
and controls dynamic range compression with a gradation curve to map the 
processed image into the appropriate dynamic range for the display. Application 
of the gradation curve will impact both the apparent detail contrast and dynamic 
range of the displayed image. Hence, both the modifying functions of the multi- 

30 resolution processing and the shape of the gradation curve affect the detail 
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contrast in the image. This requires that both be adjusted when setting the detail 
contrast in the image. 

U.S. Patent No. 5,978,518, inventors Oliyide et al., entitled "Image 
Enhancement in Digital Image Processing", issued Nov. 2, 1999 and U.S. Patent 
5 No. 6,069,979 (continuation-in-part of U.S. Patent No. 5,978,518), inventor 
VanMetter, entitled "Method for Compressing the Dynamic Range of Digital 
Projection Radiographic Images", issued May 30, 2000, describe a multi- 
resolution method for performing dynamic range modification and high-frequency 
enhancement (including detail contrast). The methods include a tone scale look- 

10 up-table that is used to map the image for display rendering. A tone scale look- 
up-table impacts the dynamic range and contrast of detail in an image. Hence, in 
this method, the dynamic range and detail contrast of the image depends on both 
the settings of the frequency modification and the parameters of the tone scale 
look-up table. It is desirable, instead to have a single set of parameters that 

1 5 control these attributes independently. 

U.S. Patent No. 6,072,913, inventor M. Yamada, entitled "Image 
Processing Method and Apparatus", issued June 6, 2000, describes a multi- 
resolution method for enhancing frequencies with dynamic range compression. 
The described invention requires the definition of many functions to control the 

20 performance of the algorithm. It does not disclose a set of parameters that directly 
and independently control all of the fundamental attributes of image quality. 

Thus, there is a need for a method that can be applied to an image, 
the parameters of which provide direct and independent control of the above 
stated fundamental attributes of image quality. 

25 

SUMMARY OF THE INVENTION 

According to the present invention, there is provided a solution to 
the problems and fulfillment of the needs discussed above by means of a method 
of digital image enhancement, especially image enhancement of medical 
30 diagnostic (radiographic) digital images. 
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According to a feature of the present invention, there is provided a 
method for enhancing an image: providing a digital image; decomposing the 
image into a multi-frequency band representation including a low frequency band 
image and multiple different high frequency band images; multiplying each of the 
5 high frequency band images with a gain factor; summing together the unmodified 
low frequency band image and the modified high-frequency band images to 
produce a reconstructed digital image: and mapping the reconstructed image 
through a tone scale look-up-table to map it to optical densities. 

ADVANTAGEOUS EFFECT OF THE INVENTION 

The invention has the following advantages. 

1 . A method for independently controlling fundamental 
attributes of image quality of a digital image. 

2. A method for controlling detail contrast, sharpness, and 
noise appearance as a function of log exposure of a digital image is provided. 

3. A method for density-dependent control of detail contrast, 
sharpness and noise appearance of a digital image is provided. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 Fig. 1 is a flow diagram of an embodiment of the present invention. 

Fig. 2 is a graphical view showing an example of tone scale 
function and parameters useful in explaining the present invention. 

Fig. 3 is a diagrammatic view illustrating the brightness control 
used in the present invention. 
25 Fig. 4 is a diagrammatic view illustrating the effect of increasing 

the slope of the tone scale useful in explaining the present invention. 

Fig. 5 is a diagrammatic view illustrating the independent control 
of detail contrast vs. dynamic range in the present invention. 

Figs. 6a-6c are graphical views of functions to control detail, 
30 contrast, sharpness, and fine detail, respectively according to the present 
invention. 



10 
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Figs. 7a-7c are graphical views of log exposure dependent 
functions to control detail contrast, sharpness, and fine detail, according to the 
present invention respectively. 

Fig. 8 is a graphical view showing the mapping of density 
5 breakpoints to log exposure break points according to the present invention. 

Fig. 9 is a diagrammatic view of density dependent control of the 
detail contrast according to the present invention. 

Fig. 10 is a flow diagram of another embodiment of the present 

invention. 

10 Fig. 1 1 is a block diagram of a digital image enhancement system 

for carrying out the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

Referring now to Fig. 1 , there is shown a block diagram of the 

15 present invention. A digital image in which code value is linearly related to log 
exposure is captured with an image acquisition unit 10. Unit 10 can be for 
example, a medical image acquisition unit such as, a diagnostic image unit (MRI, 
CT, PET, US, etc.), a computed radiography or direct digital radiography unit, an 
x-ray film digitizer, or the like. Any other digital image acquisition unit can also 

20 be used). The present invention processes the log exposure code value data, as 
shown in Fig. 1 , accordingly, the digital image data is split into four frequency 
bands B 0 (iJ), Bj(ij), B2(iJ), and Bs(iJ). The log exposure code value data I(iJ) of 
the input digital input digital image is first processed by three different low-pass 
filter operators 20, 30, 40. Each operator uses a square-wave filter. It will be 

25 evident to those skilled in the art that other low-pass filter shapes such as a 
triangle-filter can be used. The filter kernel sizes are chosen to manipulate 
different sized features (or different frequency ranges) in the image. The first 
low-pass filter operator 20 uses kernel JSTo- This operator uses the largest kernel 
size and it passes only the lowest frequencies. It generates the low-pass image 

30 Io(ij)- The next low-pass operator 30 uses kernel Ku This operator uses the next 
largest kernel size and it passes the low to mid frequencies. It generates the low- 
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pass image Ij(ij). The final low-pass operator 40 uses kernel K 2 . This operator 
uses the smallest kernel size and it passes all frequencies except for the very 
highest. It generates the low-pass image hfij)- The low-pass images are used to 
generate the frequency bands. The lowest frequency band image is Bo(iJ) and is 
5 equal to Io(iJ). This band represents large-sized features in the image (and 

contains the lowest frequencies). Manipulation of this band produces a change in 
dynamic range or latitude. The next frequency band Bi(ij) = Ij(ij) - IofiJ) is 
generated by subtractor 60. This band contains the low-to-mid frequencies and 
represents mid-sized features in the image. Manipulation of this band produces a 

1 0 contrast effect without affecting the overall dynamic range. The next frequency 
band B2(ij) = htfj) -Ii(U) is generated by subtractor 70. This band contains the 
mid-to-high frequencies and represents the small-sized features in the image. 
Manipulation of this band produces a sharpness or blurring effect of the small- 
sized features. The next frequency band B 3 (ij) =I(ij) - h(ij) is generated by 

15 subtractor 80. This band contains the highest frequencies in the image and 

represents very fine detail in the image. Manipulation of this band produces a 
sharpness or blurring effect of very fine detail in the image. 

As shown further in Fig. 1, the higher frequency bands Bj(iJ), 
B2(ij), and B$(ij) are multiplied by respective gain terms 91, 92, 93 by respective 

20 multipliers 61, 71, 81 and summed together with the lowest frequency band image 
Bo(ij) to generated the reconstructed image using V. The reconstructed image V 
is defined as 

/'(/, j) = B Q (i, j) + g } (B Q (i, J)) * B x (/, j) + g 2 (B 0 (/, j)) * B 2 (/, j) + g 3 (B 0 (/, J)) * B 3 (/, j) 

If the bands are not manipulated, that is if gi(B 0 ( g2(B 0 ( ij))= 
25 g$(B 0 ( ij))=\.Q, the reconstructed image V is equal to the original image /, i.e., 
I\iJ) = I(iJ). 

As shown in Fig. 1 the reconstructed image /' is mapped through a 
tone scale (TS) process 90 which maps the values into the desired display range. 
The output image is represented as /"(/, j) = TS(I'(i 9 j)) . 
30 Fig. 2 shows an example of a tone scale function 100 used to map 

log exposure (x-axis) to optical density (y-axis). It is a sigmoid shape and is 



parameterized by the minimum density and maximum density, the reference 
density (or pivot point), the brightness (or shift of tone scale along the log 

exposure axis relative to the pivot point), and the slope, , about the reference 

dx 

density. The reference density maps the log exposure value that will remain 
5 invariant when the slope of the tone scale is varied. The brightness control is used 
to shift the mean density of the image. The brightness control defines the log 
exposure code value that is mapped via the tone scale to the reference density. 

Fig. 3 shows the effect of the brightness control (Image A to Image 
C). Increasing the brightness control lowers the mean density of the image (Image 

10 C is lower in mean density than Image A) (or increases the brightness of the 
image), while decreasing the brightness control raises the mean density of the 
image (or darkens the image). In a preferred embodiment, the tone scale function 
is further mapped to a standard gray scale display function (as specified by the 
Medical Image Standard format DICOM PART 10) for monochrome image 

15 presentation. This facilitates the same gray scale rendering between softcopy 
(video monitor) systems and between softcopy and hardcopy (film, paper) 
systems. 

The amount of contrast in which the mid size detail in an image is 
rendered is defined as the detail contrast. In traditional screen film radiographic 

20 imaging systems, there is an inherent tradeoff between detail contrast and latitude 
(or range of log exposure that are rendered to unique densities). In these systems, 
if the contrast is increased, the latitude is decreased (or image dynamic range is 
increased). Likewise, if the latitude is increased, (or image dynamic range 
decreased) the detail contrast is decreased. The same is true if the image dynamic 

25 range and detail-contrast are controlled via a digital enhancement process (e.g., 
multi-resolution decomposition) in which the output of the process is mapped 
through a tone scale to render it to the display. Inherently, the tone scale also 
modifies both the dynamic range and detail contrast of an image. As shown in 
Fig. 4, from image A to image C as the slope of the tone scale increases both the 

30 image dynamic range and detail-contrast are increased and the latitude is 

decreased. Unless the digital enhancement process accounts for the tone scale 
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effects both the dynamic range and detail-contrast of the image will be affected 
upon display. 

In the embodiment of the present invention, shown in Fig. 1 , the 
lowest frequency band of the multi -resolution decomposition is B 0 (iJ). If this 
5 band is manipulated it affects the dynamic range of the image. As previously 

discussed, the slope of the tone scale, , also affects the dynamic range. 

dx 

According to the present invention, the lowest band is not manipulated for 
dynamic range control; instead, the slope of the tone scale is used. The dynamic 
range control is defined as the slope of the tone scale, dynamicRange = dD I dx at 

10 the reference density. Increasing the dynamic range reduces the latitude while 
decreasing the dynamic range increases the latitude. 

Manipulation of band Bj(iJ ) affects the detail contrast of the 
image. According to the present invention, band Bj(ij ) is multiplied by the gain 
term gi(B 0 (iJ)) at multiplier 61. As previously discussed, the slope of the tone 

1 5 scale also affects the detail contrast of the image. To account for the fact that both 

the slope of the tone scale and the gain term gi(B 0 (iJ)) affect the detail 

dx 

contrast, a detail contrast control is defined as a 
detailContrast(B 0 (i, j) ) = gj (B, (i, j) )*(dD/dx). 

By defining the dynamic range and detail contrast parameters in 
20 this way independent control of these attributes can be achieved. This is 

illustrated by Fig. 5. It shows a 3x3 matrix of images (i.e., IA, IB, IC, IIA, IIB, 
IIC, IIIAy IIIB, IIIC) processed using these controls. By moving from left to 
right (Image A and Image C) across a row (I, II, III) in the matrix, the dynamic 
range increases as the value of the dynamic range control is increased, while the 
25 detail contrast is held constant at a particular level. By moving from the bottom to 
top of the matrix along a column (A, B, C) the detail contrast increases as the 
value of the detail contrast control is increased while the dynamic range is held 
constant at a particular level. The image 200 (IA) in the bottom left corner has 
wide latitude and low detail contrast, while the image 201 (IIIC) in the upper right 
30 corner has narrow latitude and high detail contrast. 
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The controls for the sharpness of small features, sharpness(Bo(iJ))> 
and very fine features, flneDetail(B 0 (i j)), are defined as follows 

ft('ifcj)) 

fineDetail(B 0 ft j) ) = * 3 ffiH" 

g 2 (#oftj)) 

5 To achieve the independent control of sharpness/blurring, it is 

found that controlling the gain ratio between the frequency bands is important. 
The sharpness(Bo(ij)) parameter provides control of the sharpness of small detail 
in an image, and the fineDetail(B 0 (iJ)) parameter provides an extra level of 
control over the sharpness of very fine detail in an image. 
10 The gain terms that are used to manipulated the frequency bands, 

Bu B 2 , and B 3 are derived from detail contrast, sharpness and fine detail controls 
respectively as follows 

g x (J? 0 ft j) ) = detailContrast(B Q ft \))/dD I dx 
15 g 2 (B,(i,j)) = sharpness(B 0 (Ui))^ gl (B 0 (Uj)) 

g, (Bo ft j) ) = fineDetail(B, ft j) ) * g 2 (B 0 ft j) ) 
The present invention enables exposure dependent control of the 
detailContrast(Bo(ij)), sharpness(B 0 (iJ)),and JineDetail(Bo(iJ)). The low band 
image B 0 (iJ) represents the average log exposure and serves as a map of exposure 

20 regions in an image. To achieve exposure dependent control of the image 

enhancement, the parameters are defined as a function of that band. Areas of the 
low band image where the code values are low correspond to the low exposure 
regions in the image (that is the low density or bright areas in the image). Areas 
where the code values are high correspond to high exposure regions in the image 

25 (that is the high density or dark areas in the image). This information can be used 
to provide exposure dependent control of detail contrast, sharpness of edges, and 
sharpness of fine detail and to improve the enhancement of the image. 

Figs. 6a-6c show examples of a functional form of the detail 
contrast (Fig. 6a), sharpness (Fig. 6b), and fine detail (Fig. 6c) controls, 

30 respectively. In this example the function is constant across all exposure regions 
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(exposure independent). When a control is set to a value of 1 .0 there is no 
enhancement of features in the image. While a control value greater than 1 .0 
result in the enhancement of features in the image. Both the sharpness and the 
fine detail controls can have values less than 1 .0 (but not negative). If these 
5 controls are set to a value less than 1 .0, the corresponding features are blurred (the 
gain of these bands are decreased relative to the lower bands). Setting the 
JlneDetail(B 0 (ij)) control to a value less than 1 .0 is used to reduce the appearance 
of high frequency noise in an image. 

Figs. 7a-7c shows an example of an exposure dependent functional 

10 form that can be used for detail contrast, sharpness and fine detail controls, 
respectively. The functional form is a piecewise linear curve in which two 
breakpoints are used. The function is parameterized by setting a left and right log 
exposure break point, and a left function value (for B 0 (iJ) <= left break point) and 
a right function value (for B 0 (ij) >= right break point). The functional form 

1 5 between the breakpoints is linear. 

As shown in Fig. 8, breakpoints can be first assigned in density and 
then mapped to log exposure breakpoints via the tone scale curve. The high- 
density breakpoint is mapped to the right log exposure breakpoint and the low- 
density breakpoint is mapped to the left log exposure breakpoint. Once the 

20 brightness and dynamic range of the image are set, via the tone scale, then the log 
exposure breakpoint are defined. Setting the breakpoint in density provides a 
mechanism to achieve a preferred rendering of density regions. Often they 
correspond directly to important anatomical structures in a radiographic image of 
an individual. For example, on a chest image the lung field is a higher density 

25 (darker) region and can be rendered differently from the diaphragm, which is in a 
lower density (brighter) region. As illustrated in Fig. 9, a user or image analysis 
algorithm can independently control the detail contrast (as well as sharpness of 
small and fine detail (or noise suppression)) in low and high-density areas of an 
image. For example the detail contrast can be set higher in the lung field areas 

30 300 and lower in the lower density diaphragm area 301 . 
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It will be evident to those skilled in the art that there are other 
useful functional forms that can be generated for the detailContrast(Bo(ij)) 9 
sharpness(Bo(iJ)),and fineDetail(Bo(iJ)) controls. 

Another embodiment of the present invention is shown in Fig. 10. 
In this embodiment, the original image / from Image acquisition unit 10 and 
the low-pass images (and not the frequency band images) (from low pass filter 
operators 20, 30, 40 are directly manipulated at multipliers 400, 401, 402, and 
403, and the results are summed together at adders 501, 502, and 503 to generated 
the reconstructed image /'(i,y) which is expressed as 

I\U j) = S 3 ( 'o (h j) j) + (1 - (/ 0 (i, j) ))/ 0 (U j) + 

(g i (/ 0 (*> j) ) " gi (/o ft J) ))/, & J) + (S 2 (/ 0 ft J) ) - #3 (h ft j) ))h ft i) 

This embodiment produces the same result as the method shown in 
FIG. 1 but does not require the calculation of the frequency band images and may 
be advantageous in some implementations. 

It is a preferred embodiment of this invention that the brightness 
and dynamic range parameters and the parameters for the detailContrast(Bo(ij)), 
sharpness(Bo(iJ)),and fineDetail(Bo(iJ)) functions be either entered directly at a 
user interface to the algorithm or automatically determined via an image analysis 
method. 

Referring now to Fig. 11, there is shown a digital computer 600 for 
carrying out the present invention. As shown, digital computer 600 includes 
memory 601 for storing digital images, application programs, operating systems, 
etc. Memory 601 can include mass memory (such as hard magnetic disc or CD 
ROM), and fast access memory (such as RAM). Computer 600 also includes 
input devices 602 (such as keyboard, mouse, touch screen) display 603 (CRT, Flat 
Panel Display), central processing unit 604, output device 605 (thermal printer, 
laser printer, etc.). Components 601, 602, 603, 604, and 605 are connected 
together by control/data bus 606. Computer 600 can include a transportable 
storage medium drive 607 for reading from and/or writing to transportable storage 
media 608, such as DVD or CD. 



The invention has been described in detail with particular reference 
to certain preferred embodiments thereof, but it will be understood that variations 
and modifications can be effected within the spirit and scope of the invention. 
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PARTS LIST 



1 0 image acquisition unit 

20 low-pass operator using filter Ko 

30 low-pass operator using filter Kj 

40 low-pass operator using filter K2 

60.70.80 subtracter 

61.71.81 multiplier 

62.72.82 adder 

90 tone scale process 

91,92,93 gain terms 

1 00 example tone scale function 

200 low detail contrast and wide latitude image example 

201 high detail contrast and narrow latitude image example 
300 high-density lung field area 

3 0 1 low-density diaphragm area 

400-403 multiplier 

501-503 adder 

600 computer 

601 memory 

602 input device 

603 display 

604 processing unit 

605 output device 

606 control/data bus 

607 transportable storage medium 

608 storage media 



